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ABSTRACT 

The globular cluster HP 1 is projected on the bulge, very close to the Galactic 
center. The Multi-Conjugate Adaptive Optics (MCAO) Demonstrator (MAD) 
at the Very Large Telescope (VLT) allowed to acquire high resolution deep im- 
ages that, combined with first epoch New Technology Telescope (NTT) data, 
enabled to derive accurate proper motions. The cluster and bulge field stellar 
contents were disentangled by means of this process, and produced unprecedented 
definition in the color-magnitude diagrams for this cluster. The metallicity of 
[Fe/H]f«-1.0 from previous spectroscopic analysis is confirmed, which together 
with an extended blue horizontal branch, imply an age older than the halo aver- 
age. Orbit reconstruction results suggest that HP 1 is spatially confined within 
the bulge. 

Subject headings: globular clusters: general — globular clusters: individual: HP 1 
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Introduction 



A deeper understanding of the globular cluster population in the Galactic bulge is 
becoming possible, thanks to high resolution spectroscopy and deep photometry, using 
8-10m telescopes. An interesting class concerns moderately metal-poor globular clusters 
([Fe/H]«-1.0), showing a blue horizontal bra nch (BHB). A dozen of these objects are found 



projected on the bulge ( jBarbuy et al. 



20091 ). and they might represent the oldest globular 



clusters formed in the Galaxy. 
According to iGao et al. 



( 120101 ) . the first generation of massive, fast-evolving stars, 
formed at redshifts as high as z^35. Second generation low mass stars would be found 
primarily in the inne r parts of the Galaxy today, as well as inside satellite galaxies. 



Nakasato fc Nomoto 



(120031 ) suggest that the metal-poor component of the Galactic bulge 
should have formed through a subgalactic clump merger process in the proto-Galaxy, 
where star formation would be induced and chemical enrichment by supernovae type II 
occurred. The metal-rich component instead would have formed gradually in the inner disk. 
Therefore, metal-poor inner bulge globular clusters might be relics of an early generation of 
long-lived stars formed in the proto-Galaxy. 

Another aspect of the interest of inner bulge studies, comes from evidences that stellar 
populations in the Galactic bulge are similar to those in spiral bulges and elliptical galaxies, 



(Bica 


1988; 


Rich 


1988) 



19881 ) . In the past, the detailed study of the bulge globular clusters was 
hampered by high reddening and crowding. With the improvement of instrumentation, it is 
now possible to derive accurate proper motions, and to apply membership cleaning. Most 
previous efforts in the direction of proper motio n studies were carrie d out using Hubble 



Space Telescope - H S T dat a by our group (e.g. ( 



Feltzing &: Johnson 



Zoccali et al. 



feOOJ ) for NGC 6528, and kuijken fc Rich 



2001)) for NGC 6553, or 



( 120021 ) for field stars, 
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among others. 



In a systematic study of bulge globular clusters (e.g. iBarbuy et al. I ( 119981 . 120091 ). we 
identified a sample of moderately metal-poor clusters ([Fe/H]~-1.0) concentrated close 
to the Galactic center. HP 1 has a relatively low reddening, and appeared to us as a 
suitable target to be explored in detail in terms of proper motions. Its coordinates are 
a = 17 h 31 m 05.2 s , 5 = -29°58'54" (J2000), and projected at only 3.33° from the Galactic 
center (l=-2°.5748, b=2°.1151). 

In 2008, the European Southern Observatory (ESO) announced a publ ic call for Science 



Marchetti et al. 



Verific ation of the Multi-Conjugate Adaptive Optics Demonstrator (MAD 
(120071 )). installed at the ESO UT3-Melipal telescope. Several star clusters were observed 
with the MAD facility, deliv ering high quality data for the open clusters Trapezium, 



FSR 1415 and Trumpler 15 (IBouy et al. 



2008 



the globular clus ters Terzan 5 anc 



and 30 Doradus (ICampbell et al. 



Momanv et al 



NGC 3201 (IFerraro et al. 



2008 



20091 : 



Sana et al. 



Bono et al. 



2010) 



20101 ) 



20101 ). The major advantage of MAD is that it allows 
correcting for atmospheric turbulence over a wide 2' diameter field of view, and as such, 
constitutes a pathfinder experiment for future facilities at the European Extremely Large 
Telescope (E-ELT). Wide field adaptive optics with large telescopes opens a new frontier 
in determining accurate parameters for most globular clusters, that remain essentially 
unstudied because of high reddening, crowding (cluster and/or field), and large distances. 

In Sect. 2 the observations and data reductions are described. In Sect. 3 the 
HP 1 proper motions are derived. The impact of the high quality proper motion cleaned 
Color-Magnitude Diagrams (CMDs) on cluster properties are examined in Sect. 4. The 
cluster orbit in the Galaxy is reconstructed in Sect. 5. Finally, conclusions are given in 
Sect. 6. 
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2. Observations and Data Analysis 

MAD was developed by the ESO Adaptive Optics Department to be used as a visitor 
instrument at VLT-Melipal in view of an application to the E-ELT. It was installed at 
the visitor Nasmyth focus and the concept of multiple reference stars for layer oriented 
adaptive optics corrections was introduced. This allows a much wider and more uniform 
corrected field of view, providing larger average Strehl ratios and making the system a 
powerful diffraction limit imager. This is particularly important in crowded fields where 
photometric accuracy is needed^. Following our successful use of MAD in the first Science 
Demonstration (Momany et al. 2008), we were granted time to observe HP 1 in the second 
Science Demonstration (p). 

HP 1 was observed on August 15 th , 2008. Table [6] displays the log of the J and K s 
(for brevity we use K) observations. Clearly, the seeing in K was excellent being almost 
half that in J (0'.'23 vs 0'.'38). The MAD infrared scientific imaging camera is based on a 
2048 x 2048 pixel HAWAII-2 infrared detector with a pixel scale of 0"028. In total, 25 min. 
of scientific exposures were dedicated to each filter, and subdivided into 5 dithered images. 

The images were dark and sky-subtracted and then flat-fielded following the standard 
near infrared recipes (e.g. Momany et al. 2003), within iraf environment. The typical 
dithering pattern of a MAD observation allows a field of view of ~ 2' in diameter. Within 
this field of view, three reference bright stars (R magnitudes that ranged between 12.5 
and 13.8 according to their UCAC2 magnitude system) were selected to ensure the optics 
correction. However, one of these proved to be a blend of two stars, which did not allow a 
full optical correction of the field. 



1 ht t p : / / www . eso . org/pro j ect s/ aot / mad / 
2 http: / / www.eso.org/ sci / activities / vltsv / mad / 
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Figure [T] shows the mosaic of all 10 J and K images as constructed by the daophot- 
/montage2 task. The superb VLT/MAD resolution is illustrated when compared to that 
seen in the 2MASS K image of HP 1. For the proper motion purposes, the MAD data set 
represents our second epoch data. The displacement of the HP 1 stellar content was derived 
by comparing the position o f the stars in this data set with respect to that of SUSI@NTT 
obtained on May 16' ft , 1994 ( Ortolani et al. Ill997 ). The epoch separation is 14.25 years. It 
is worth emphasizing that the seeing of 0.4 5 " for the V image of HP 1 is one of the best 



obtained at the NTT (lOrtolani et al. 



1994 



19971 ) 



The photometric reductions of the two data sets (NTT and MAD) were carried out 
separately. We found about 3100 stars in common between the MAD and NTT data, which 
were used for the proper motion analyses. 



2.1. Photometric Reduction and Calibration 



The stellar photometry and astro metry were obtained by point spread function (PSF) 



fitting using daophot ii/allframe ( jStetson et al. 



19941 ). Once the FIND and PHOT 



tasks were performed and the stellar-like sources were detected, we searched for isolated 
stars to build the PSF, for each single image. The final PSF was generated with a PENNY 
function that had a quadratic dependence on frame position. ALLFRAME combines PSF 
photometry carried out on the individual images and allows the creation of a master list by 
combining images from different filters. Thereby this pushes the detection limit to fainter 
magnitudes and provides a better determination of the stellar magnitudes (given that 5 
measurements were used for each detected star). 

Our observing strategy employed the same exposure time for all images, and no bright 
red giant stars were saturated. When producing the photometric catalog in one filter (and 
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since only the central part of the field of view had multiple measurements of any star), 
stars appearing in any single image were considered to be real. Later, when producing the 
final JK color catalog, only those appearing in both filters were recorded (this way we 
removed essentially all detections due to cosmic rays and other spurious detections). The 
photometric catalog was finally transformed into coordinates with astrometric precision by 
using 12 UCAC2 reference starjf] with R < 16.2. 

Photometric calibration of the J and K data has been made possible by direct 
comparison of the brightest MAD non-saturated (8.0 < K < 13.0) stars with their 2MASS 
counterpart photometry. 

From these stars we estimate a mean offset of A Jj@mad-j@2MASS — —3.157 ± 0.120 
and A K k @mad-k@2MAss = -1-395 ± 0.134. 

Photometric errors and completeness were estimated from artificial star experiments 
previously applied to similar MAD data (Momany et al. 2008). 

The images with added artificial stars were re-processed in the same manner as 
the original images. The results for photometric completeness showed that we reach a 
photometric completeness of ~ 75%, 50% and 10% around K ~ 17.0, 17.5 and ~ 18.0, 
respectively. 

3. Derivation of the proper motions 

The proper motion of the HP 1 stellar content was derived by estimating the 
displacement in the (x, y) instrumental coordinates between MAD (second epoch data) and 
the NTT (first epoch) data. Since this measurement was made with respect to reference 



3 http:/ /ad. usno.navy.mil/ucac/ 
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stars that are cluster members, the motion zero point is the centroid motion of the cluster. 
The small MAD field of view is fully sampled by the wider NTT coverage, and thus 
essentially all MAD entries had NTT counterparts for proper motion determination. In this 
regard, we note that the photometric completeness of the MAD data set is less than that 
of the optical V, I NT T data, that reached at least 2 magnitudes below the cluster turnoff 



( jOrtolani et al. 



19971 ). 



The daophot/daomaster task was used to match the photometric MAD and NTT 
catalogs, using their respective (x, y) instrumental coordinates. This task employed a 
cubic transformation function, which, with a matching radius of 0.5 pixels or 0'.'015, easily 
identified reference stars among the cluster's stars (having similar proper motions). As a 
consequence, the stars that matched between the 2 catalogs were basically only HP 1 stars. 
In a separate procedure, the J2000 MAD and NTT catalogs transformed into astrometric 
coordinates, were applied to a matching procedure, using the iraf/tmatch task. This 
first/second epoch merged file included also the (x,y) instrumental coordinates of each 
catalog. Thus, applying the cubic transformation to both coordinate systems yielded the 
displacement with respect to the centroid motion of HP 1. An extraction within a 1.5 pixel 
radius around (0, 0) in pixel displacement showed to contain most, and essentially only 
HP 1 stars, and we will use this selection for the rest of the paper. 

On the other hand, stars outside this radius (i.e. representing the bulge populations) 
are more dispersed, and required a careful analysis. In a first attempt we used only the 
stars with Ay < -1 (in order to avoid cluster stars) to get the baricentric position of the 
field bulge population in x, and stars with A x > 1 to measure the field bulge population 
in y. After a number of tests with different selections, in a second attempt we made a two 
gaussian component fit along a 0.2 pixels wide strip connecting the center of the cluster 
distribution and the field distribution. For the present analysis we adopted for the proper 
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motion of the cluster, relatively to the field, the mean of the two determinations (A x, A 
y) = (-2.1, 1.96). From the gaussian fit we derived the width of the distributions of the 
cluster and field stars, resulting respectively <r=0.583 and 2.565 pixels, corresponding to 
15 and 41 mas respectively. This is the error of the single star position. The statistical 
error of the mean is obtained by dividing the single star error by the square-root of the 
number of stars used (about 1700 in the cluster and 1500 in the bulge) for each of the 
two measurements. The error on the mean position of the field obviously dominates and 
corresponds to about 0.05 pixels, or 1.3 mas (about 0.1 mas/year). This is indeed a very 
small error. We performed a number of tests fitting the cluster stars and the field stars in 
different conditions checking for additional uncertainties. A further quantitative evaluation 
of the fitting uncertainty can be obtained measuring the different distances between the 
centers of the HP1 proper motions cloud to that of the field, propagated from the error on 
the angle of the line connecting the two groups of stars. The angle of the vector joining the 
center of the HP1 proper motion cloud to that of each star of the bulge was computed, and 
a Gaussian fit was performed which returned a 1-sigma dispersion of 14°. This uncertainty 
propagates on the distances of the two groups by about 0.39 mas/year, dominating over the 
other discussed sources. 



3.1. Astrometric errors 



T he astrometric errors a re the combination of different random and systematic 



errors ( 1 Anderson et al. 



20061 ) . They concluded that in the case of relative ground-based 
astrometry in a small rich field, the main error sources are random centering errors due 
to noise and blends, and random-systematic errors due to field distorsions, and finally by 



systematic errors due to chromatic effects. 
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i) Centering errors: Considering that we used two very different data sets, a first epoch 
based on direct CCD photometry, and a second one, on the corrected infrared adaptive 
optics (with a much wider scale and better resolution), we can assume that the astrometric 
random errors are dom inated by the first epoch classical CCD photometry. Following 



Anderson et al. 



(120061 ). we performed a test with 3 consecutive NTT images of Baade's 
Window taken in the same night. They were obtained with equal (4 min.) exposure times, 
and nearly at the same pixel positions. The seeing was stable and these images were 
obtained about 1 hour later than the 10 min. cluster V exposure. 

This Baade's Window field was chosen because the stars are uniformly distributed 
across the image, and the density of st ars is similar to the HP 1 field. These images 



were analysed in 



Ortolani et al. 



( 119951 ) . The position of the stars in common were 
compared. The photometric and astrometric errors are shown in Figs. HJ The astrometric 
centering error, averaged over the whole dynamic range, is 0.1 pixels, as indicated in Fig. 
2\ This corresponds to 13 mas, which is comparable with the 7 mas pointing error by 



Anderson et al. 



(120061 ). The value of 13 mas amounts to 0.91 mas/yr in an interval of 
14.25 yr. Taking into account the number of about 3000 stars (half field and half cluster) 
used for the proper motion derivation of HP 1, this leads to a minor error of 0.022 mas. 

ii) Field distortion errors: The field distortion analysis is usually performed employing 
shifted images of a field with uniformly distributed stars. HP 1 multiple shifted images 
are not available in our first epoch run. To estimate the distortions, we compared the 
proper motion of cluster stars in HP 1 NTT images as a function of the distance from the 
optical center. The distribution did not show any relevant trend, and the statistical analysis 
indicates that there is no significant shift above 0.01 pixel across the field. This corresponds 
to 0.19 mas/yr. 

Hi) Chromatic errors: Chromatic errors are due to the different refraction, both 
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atmospheric and instrumental, on stars with different colors. The refraction dependence on 
wavelength (inverse quadratic to first approximation) makes this effect more pronounced in 
V, I than in the near infrared. 

Ideally, a chromatic experiment should include a measurement of the displacement 
of stars with different colors, at different airmasses. However, we do not have spec i fic 



Anderson et al. 



(120061 ) 



observations for this test. Thus, we followed the procedure given by 
where they measured the displacement of stars as a function of their colors. We took the 
resulting pixel displacement of the cluster stars, and separately checked the shift variations 
with the color (V-I). Such plot does not show an evident dependence on color. This is 
expected since our NTT observations were taken very close to the zenith ( air mass ~ 1.04), 
and the infrared MAD data have a negligible effect. 

In order to further quantify any chromatic systematic effect, we subdivided the sample 
of cluster stars into a redder (V-I > 2.2) and a bluer (V-I < 2.2) groups, of 106 and 171 
stars respectively. The displacement between the two groups resulted to be 0.031±0.03 
pixels, corresponding to 0.9 mas. In 14.25 yr, this gives 0.06 mas/yr. 

iv) Total errors: 

By quadratically adding the errors on centering, distortion and chromatic effects of 
respectively 0.022 mas/yr, 0.19 mas/yr, and 0.06 mas/yr, a total contribution of these 
errors of 0.2 mas/yr is obtained. 

The estimated error of 0.39 mas/yr in the proper motion value indicates that the effect 
due to the mutual field and cluster stars contamination dominates over the astrometric 
pointing, distortion and chromatic errors. 
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4. The proper motion cleaned Color-Magnitude Diagram of HP 1 



Figure [3] presents the MAD K vs. J — K CMDs showing the proper motion 
decontamination process in three panels: for the whole field, cluster proper motion 
decontaminated, and remaining field only. In the top panels are shown the displacements 
of stars between the NTT and MAD images, plotted in the MAD pixel scale. 

The decontaminated cluster CMD provides fundamental parameters to study its 
properties, in particular metallicity and age. Figure H] shows a J vs. V — K proper motion 
decontaminated CMD of HP 1. A distance modulus of (m-M)j = 15.3 an d a reddening of 



20081 ) with metallicity 



E(V-K) = 3.3 were applied, and a Padova isochrone (IMarigo et al. 
Z = 0.002 and age of 13.7 Gyr is overplotted. The fit c onfirms the cluster m etallicity of 
[Fe/H]= —1.0, found from high resolution spectroscopy (IBarbuy et al. 1 120061 ) . We note 
that the turnoff is not as well matched as the giant branch, due to two well-known main 
reasons: the bluer turnoff in isochrones relative to observations is possibly connected with 
color transformations, and on the other hand, systematic bias in the photometric errors 
give brighter magnitudes close to the limit of the photometry. Similarly, Fig. [5] displays 
the K vs. V — K HP 1 di agram as compared with the NGC 6752 ([Fe/H]= —1.42) mean 



locus ( jValenti et al. 



20041 ) . The red giant branches of M 30, M 107, 47 Tuc and NGC 6441 
-1.91,-0.87,-0.70, and —0.68, respectively, are al so overplotted. 



of metallicities [Fe/H] = 
where the metallicity values, mean loci and the red giant fiducials are from 



Valenti et al. 



( 120041 ) . The HP 1 bright red giants clearly overlap the M 107 fiducial (reflecting their 
similar metallicity) and are redder with respect to the NGC 6752 bright giants. On the 
other hand, the cluster old age is reflected by the presence of a well defined and extended 
BHB (very similar to NG C 6752). Five R R Lyra e cand idates appear in the RR Lyrae gap, 



at 4.2 < (V — K) < 5.0. flTerzan 



1964a 



b, 



1965 



19661 ) reported 15 variable stars in HP 1, 



but none has been identified as RR Lyrae. The Horizontal Branch (HB) morphology is 
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sensitive mainly t o metallicity and age. T 



1994) 



param e ter effect (ISandage &: Wildey 



Rev et al. 



Potter et al. 



le age effect is related to the so-call ed second 
1 demonstrated in m odels by 



19671 ). as we 



(l200lh and in observations by 



Lee et al. 



Potter et al. 



2010|). In particular, 



(120101 ) analysed the HB morphology from ACS/HST observations, based on 
the difference between the average HB V — I color and the subgiant branch (SGB) color 
A V 

~ Ihb B i an d concluded that age dominates the second parameter. This indicator is 
very sensitive to the cluster age, and more so around the metallicity Z = 0.002. 

For HP 1 the mean V — I color difference between the SGB and the HB is 
AV- lfj G B B = 0.75 ± 0.01. A comparison with 



Potter et al. 



(j2010f )'s sample shows that 



HP 1 has a very blue HB for its metallicity. We selected 5 clusters with the same HB 
morphological index, and found an average metallicity of [Fe/H]= —1.9 ± 0.36, and another 
group of 5 clusters, with a comparable metallicity to HP 1. This second group has an 
average A V — I = 0.46 ± 0.27 which is considerably smaller than in HP 1, consistent with 
HP 1 having a much bluer HB. Both cluster groups have a mean age of 12.7 ± 0.4 Gyr. 
From Fig. 17 of Potter et al., we get an age difference of about 1 Gyr older for HP 1 relative 
to their sample of halo clusters with ~12.7 Gyr, resulting for HP 1 an age of ~ 13.7 Gyr. 
Therefore HP 1 appears to be among the oldest globular clusters in the Galaxy. 

For the distance determination, there are basically two methods: a) the relative 
distance between the cluster and bulk of the bulge field, and b) based on the absolute 
distance, which requires reddening values. For the first of these methods, we rely on the 
difference between the HP 1 horizontal branch at the RR Lyrae level, at V = 18.7, and that 
of the bulge field at V = 19.35. Thus the cluster is A V = 0.66 brighter than the field , and 



taking into account metallicity effects on the HB luminosity ( ( jBuonanno et al. 1119891 ). we 
obtain A V — 0.35 ± 0.14. This implies that the cluster is 1.2±0.4 kpc in the foreground 
of the bulge bulk population. The uncertainty is due to the metallicity difference of about 
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1 dex between the bulge ([Fe/H]«0.0) a nd HP 1 ([Fe/Hl^-l.O ) 



the Galactic center to be Rgc=8-0±0.6 (IMajaess et al 



2009 



Assuming the distance o f 



Vanhollebeke et al. 



20091 ) 



a distance of d = 6.8 kpc is obtained. If the orbits of stars near the super massive black 



hole near the Ga 



Gillessen et al. 



actic center method is used, a distance of Rgc=8.33±0.35 kpc is given by 



((20091). In this case the distance of HP 1 to the Sun is d Q = 7.1 kpc. This 
relative distance method is reddening independent, because the reddening of cluster and 
surrounding field is expected to be the essentially the same, due to a negligible reddening 
inside the bulge (Barbuy et al. 1998). Therefore the distance of the cluster to the Galactic 
center depends only on the assumed distance of the Galactic center. 

The second method of absolute distances requires reddening determinations. From the 
optical and infrared CMDs, and adopting the absolute-to-selective absorption R v = 3.2 



tearbuy et al. 



19981 ). we obtain a mean distance from the Sun of d = 7.3 ± 0.5 kpc. 
Within the uncertainties for the cluster and Galactic center distances, we conclude that 
HP 1 is probably the globular cluster located closest to the Galactic center. 

For simplicity we assume the distance of HP 1 from the Sun to be d = 6.8 kpc 
hereafter. 



5. Spatial motion of HP 1 in the Galaxy 



5.1. Absolute proper motion 



To compute the velocity components of HP l's motion we need its radial velocity 



and the proper motion. The heliocentric rad ial velocity v^ eL = 45 .8 ± 0.7 km s 1 was 



Barbuy et al. 



(120061 ). The proper motion 



adopted from the high resolution analysis by 
can be computed with respect to the bulge, and the bulge proper motion can then be 
subtracted. The bulge proper motion is a composition of the bulge internal kinematics, 
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and the reflected motion of the LSR. Near the position of HP 1, the bulge kinematics is 
close to that of a rotating solid body with spin axis orthogonal to the Galactic plane (e.g. 



((Zhao. 



19961 ) ). Note also that this correction is not large: HP 1 is projected so close to the 



Galactic center that the rotational velocity of the bulge is closely aligned with the line of 
sight, so its tangential component is very small. Because both motions are parallel to the 
Galactic plane, it is convenient to work in Galactic coordinates. 

As Fig. [6] shows, the HP1 motion with respect to the bulge has a vector = (—2.10- 
px, 1.96 px^ 



To estimate the bulge internal motion, we note that 



Tiede et al. 



(119991 ) obtained a 



rotational velocity of ~ 75km s l , in their fields 6/05 and 6/18 (b=-6° and field number), 



at positions I = —8°. 7 and I = 8°. 4, respectively. 
f rot depends linearly on the radius r (confirmed by 



n the solid body ap proximation 



Howard et al. 



(|2008|)), and at the 



4 Note that in §3] we were working with NTT — MAD coordinates, while here we put the 
two epochs in chronological order, thus MAD — NTT coordinates are used. The signs are 
therefore inverted. Taking into account the scale of 0'.'028 px _1 and that ARA = —Ax, the 
vector expressed in Equatorial coordinates is = (0'.'0588, 0'.'05488). The position angle of 
the b axis with respect to the DEC axis is —56°. 745, which can be used to rotate «iq obtaining 
UGai = (Al, Ab) with Al = 0T0781 ± 0'.'0056 and Ab = -0"0191 ± 0'.'0056). T he coordinate 



transformations were performed with the code SM (ILupton &: Monger 



19971 ). in particular 



using the package project, written by M. Strauss and R. Lupton. With an epoch difference 
At = 14.2519 ± 0.0027 yr, these displacements yield proper motion components relative to 
the bulge /i[ el x cos b ~ 5.5 mas yr -1 and /i£ cl ~ —1.34 mas yr -1 . As recalled above, to 
obtain absolute proper motions we subtract the reflected motion of the LSR, and the bulge 
internal motion. Assuming Vlsr — —243 km s _1 , at a distance of 8.4 kpc from the Galactic 
center (these values are explained below) we obtain (/xj, /j&)lsr — (6.102,0) mas yr -1 . 
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position of HP1 (Sect. 1) we expect v mt ~ 75km s _1 / 8 - 5 ° x 2 °.58 = 22.8 km sec" 1 
(where we assumed I oc r for small values of I). Basically the same value (22.7km s" 1 ) 
is predicted at a distance of 0.38 kpc (= r x sin 2°. 58 with r = 8.4 kpc) by the 
angular velocity of 60.0km s~ 1 kpc~ 1 adopted in our model (Sect. 15. 2\i . Assuming 
that bulge stars are at the distance of the Galactic Center, the tangential component 
of f rot is t>T — 22.8km s" 1 x sin 2°. 58/ cos 2°. 58 = 1.03 km s -1 directed toward 
increasing I. After adding vt to Vlsr, the composite motion of the bulge becomes 
(^i, A i fe)bui g e+LSR = (6.127, 0) mas yr -1 , and subtracting it from the relative motion quoted 
above yields for HP 1 ^ x cos b = —0.65 ± 0.39 mas yr -1 and fi^ = —1.34 ± 0.39 mas yr -1 . 
Using this proper motion, in the next section we calculate the cluster orbit. 



5.2. HP 1 orbit in the Galaxy 



he orbit was computed both with the axisymmetric model by 



Allen &: Santillan 



( 119911 ) and with a model including a bar. The models and integration algorithm are 



described in Ji l kova et al. (in preparation). Similar models have also been used in 



Magrini et al. 



(120 10[ l. Compared to these earlier versions, the axisymmetric model was 



rescaled to match the more recent valu es of rotation velo city, solar Galactocentric distance 



and solar velocity relative to the LSR. 



Reid et al. 



(120091 ) estimated a rotation velocity of 



254 ± 16 km s 1 , a nd a distance of 8.4 ± 0.6 kpc, using the solar motion relative to the 
LSR determined by 



Dehnen & Binnev 



recently re-examined by 



Schonrich et al. 



(1998). However, the analysis of Dehnen et al. was 



(120101 ) obtaining slightly different values - the 



component in the direct ion of solar Galactic rota tion is 7 km s 1 higher. Taking this into 



account, we rescaled the 



Allen fc Santillan 



( 1199 ll ) parameters to get va l ues o f 243 km s 1 at 



8.4 kpc, which is also consistent with the results of 



Reid fc Brunthaler 



( 120041 ) who obtained 



the solar rotation velocity from the proper motion of Sgr A* 
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The Galactic bar is modeled by a Ferrers potential of an inho mogeneous triax ial 



ellipsoid (IPfenniger 



Pichardo et al. 



(120041 ) 



19841 ) . The model parameters are adopted from 
with length of 3.14 kpc, axis ratio 10:3.75:2.56, mass of 0.98 x 10 10 M , angular velocity of 
60.0 km s~ 1 kpc~ 1 , and an initial angle with respect to the direction to the Sun of 20° (in 
the direction of Galactic rotation). For the axisymmetric background we keep the potential 
described above with decreased bulge mass by the mass of the bar. 

The initial conditions for the orbit calculations are obtained from the observational 
data characterizing the cluster: coordinates, distance to the Sun, radial velocity, and proper 
motion. To evaluate the impact of the measurement errors, we calculate a set of 1,000 
orbits with initial conditions given by sampling the distributions of observational inputs. 
We assume normal distributions for the distance to the Sun, radial velocity, and proper 
motion. The errors on radial velocity and proper motion components are given above, while 
for the distance to the Sun we assumed an error of 10%. 

The transformation of the observational data to the Cartesian c o ordin ate system 



centered on the Sun was carried out with the 



Johnson fc Soderblom 



(119871 ) algorithm. 



The velocity vector with respect to the LSR i s then obtained by cor rection for 



the solar motion with respect to the LSR from ISchonrich et al. 



J2010): (U,V,W) Q 



(11.1, 12.24, 7.25) km s _1 (right-handed system, with U in the direction to the Galactic center 
and V in the Galactic rotation direction). The final transformation to the Galactocentric 
coordinate system was made by using the solar Galactocentric distance 8.4 kpc and the 
LSR rotation velocity of 243km s _1 (see above). We obtain proper motions in equatorial 
coordinates /ira x cosDEC = 0.76 ± 0.39 mas yr^ 1 , /xdec — —1-28 ± 0.39 mas yr _1 , and 
Cartesian coordinates and velocities (x,y,z) = (1.59, 0.31, 0.25) kpc and (v x , v y , v z ) = 
(—57.47, —231.99, —34.26) km s -1 . We adopted a right-handed, Galactocentric Cartesian 
system (x to the Sun direction, z to the north galactic pole). 
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We integrate the orbits with such initial conditions b ackwards for an in terval of 3 Gyr 



19921 ). The example 



using a Bulirsch-Stoer integrator with adaptive time-step (IPress et al. 
of orbits given by the average values of the observational data is given in Fig. [7J The 
presence of the bar disturbs the orbit of this central globular cluster, causing deviations 
not found in the axisymmetric model. This can be considered as an upper limit for the 
excursions that the cluster can make inward and outward. Even so, it is clear that the 
cluster is essentially confined within the bulge. 



Running simulations for a longer time is not very mea ningful because t 



lere is evidence 



Minchev et al. 



( 120101 ) suggest 



that the bar structure is a transient feature. For example 
that the current bar might have formed only 2 Gyr ago. It is then impossible to simulate 
the orbit along the entire life of the Milky Way, but one can guess that older bars would 
have had a similar effect on the orbit of HP1. Note that we also included the spiral arms, 
but they do not change the orbit significantly. They are very weak and the orbit is too close 
to the Galactic center to be influenced by any radial migration due to bar and spiral arm 
interaction. 

We calculated orbital parameters as averaged values over individual revolutions in 
the Galactic plane for each orbit. Distributions for apogalacticon R a , vertical height of 
orbit |z m ax|, arid phase period T are shown in Fig. [HI In general the orbits do not reach 
galactic distances larger than 5 kpc and also the cluster remains close to the Galactic plane 
(kmaxl < 0.3kpc for the axisymmetric model, |z max | < 0.6kpc for the model including bar). 

For compa r ison p urposes, we also computed the HP 1 orbit using the code developed by 



Mirabel et al. 



(120011 ) . that includes the Galactic s pheroidal and the di sk potentials. It was 



recently applied to u Centauri orbital simulations (jSalerno et al. 



20091 ) . In this method we 



use essentially the same initial conditions (U ,V ,W ), as in the previous method above, 
and the simulation results agreed well with the previous method for the barless model. 
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Conclusions 



The clear definition of an extended blue horizontal branch morphology obtained from 
these high spatial resolution data, as provided by the proper motion cleaning method, 
indicates a very old age for HP 1, of ~1 Gyr older than the halo average. 

The proper motions and orbits derived indicate that HP 1 does not wade into the halo 
and is confined within the Galactic bulge. As a consequence, HP 1 can be identified as 
a representative relics of an early generation of star clusters formed in the proto-Galaxy. 
The very old globular cluster N GC 6522, also haying m oderate metallicity and BHB, is as 



well confined within the bulge (ITerndrup et al 



metal-rich bulge globular cluster NGC 6553 ( IZoccali et al. 



1998 ). As compared with the templat e 



2001 



Ortolani et al. 



19951 



HP 1 appears to have a more excentric orbit, and it is much closer to the Galactic center. 

Extensive tests of orbits within potential wells that includes a massive bar show that 
the confinement of HP 1 within the bulge is maintained, even in case of random orbits 
generated by the presence of the bar. 

The case of HP 1 with wide field multi-conjugate adaptive optics shows that such 
ground-based facilities can be used for high spatial resolution studies of crowded inner bulge 
clusters. Such data can provide a great impact for a better understanding of the globular 
cluster subsystems, and their connection with stellar populations in the Galaxy, and the 
sequence of processes involved in the formation of the Galaxy itself. 



We are grateful to Ata Sarajedini for helpful discussions on the second parameter effect. 
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the ESO Adaptive Optics group, in particular, Enrico Marchetti and Paola Amico. BB and 
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the Italian Ministero de l'Universita e della Ricerca Scientifica e Tecnologica. 
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Table 1: Log of the MAD observations obtained on August 15 th 2008. 
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Fig. 1. — MAD image of the globular cluster HP 1, obtained at the ESO-VLT at Paranal, 
Chile. Left panel is a composite color image of HP 1 from the near-infrared 2MASS survey; 
Right panel is a close-up of the MAD i^-band image covering the central If 8 x 1'.8 of HP 1. 
North is up, East to the left. 
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Fig. 2. — Photometric and astrometric errors of Baade's Window images. Upper 
panel:Photometric errors AV vs. instrumental magnitude. Dashed line shows the aver- 
age astrometric error. Lower panel: Astrometric error in NTT pixel scale vs. instrumental 
magnitude. 
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Fig. 3. — HP 1: Proper motion decontaminated K vs. J — K CMD. Upper panels: Left: 
Displacement of NTT (1st epoch) - MAD (2nd epoch), plotted in MAD pixel scale; Middle: 
Cluster sample is encircled where I, b directions are indicated by the arrows; Right: Field 
sample only. Lower panels: Left: Observed CMD; Middle: Proper motion decontaminated 
cluster CMD. The two stars indicated as red squares are the two giant stars analysed with 
high resolution spectroscopy; Right: Remaining field stars CMD. The extraction is within a 
1.5 pixel radius. 
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V-K 



Fig. 4. — HP 1: J vs. V — K CMD, built from a proper motion selected subsample of 
stars. The extraction is within a 1.5 pixel radius, in pixel displacement. Padova isochrone 
of Z = 0.002 and age of 13.7 Gyr is overplotted on the observed CMD. A distance modulus 
of (m-M)j = 15.3 and a reddening E(V-K) = 3.3 were adopted. The metallicity, distance, 
and reddening adopted for the fit are in agreement with the spectroscopic analysis (Barbuy 
et al. 2006). The present MAD photometry reaches the turn-off limit. 
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Fig. 5. — The K vs. V — K HP 1 diagram (black dots) as compared with the NGC 6752 
(green dots) catalog (Valenti et al. 2004). We assumed the fiducial values (to— M) = 13.18, 
E(B — V) = 0.04, Av = 3.1, Ak = 0.38. Also plotted are the red giant branches of 
comparison Galactic globular clusters (blue lines): M 30, M 107, 47 Tuc and NGC 6441 
([Fe/H]— 2.12, —1.04, —0.76 and —0.68 (where the catalog and the red giant fiducials are 
from Valenti et al. 2004). The HP 1 bright red giants are clearly overlapping the M 107 
fiducial (reflecting their similar metallicity) and are redder with respect to the NGC 6752 
bright giants. 
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Fig. 6. — In the reference system (Ax, Ay) = (xmad — ^ntt,?/mad — Vntt) with origin in 
HP1 the bulge has coordinates (2.10,-1.96). Therefore the HP1 motion with respect to 
the bulge has a vector = (—2.10, 1.96) (red arrow). Because RA and Ax have opposite 
directions, the vector expressed in Equatorial coordinates is «iq = (0'/0588, C'05488). The 
relative orientation of the Equatorial and Galactic coordinate systems is shown in the upper 
right corner, making clear that the proper motion of HP1 is mostly along positive I, with a 
small component along negative b. The orbit must then be confined near the Galactic plane. 
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Fig. 7. — Galactic orbit of the cluster: projections into the Galactic and meridional plane 
are plotted on the left and right, respectively. Orbit in axisymmetric model is plotted by red 
line, orbit in the model including bar by blue line. The initial conditions are given by mean 
observational input data. See text for detailed description of integrations. 
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Fig. 8. — Distributions of orbital parameters: apogalacticon R a , vertical height of orbit 
l^maxl; an d phase period T. Distributions given by the axisymmetric model are plotted by 



red line and in the model with bar by blue line (same as in Fig. [7j). 



